The role of initial state reconstruction in short and long time deviations 

from exponential decay. 
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We consider the role of the reconstruction of the initial state in the deviation from exponential 
decay at short and long times. The long time decay can be attributed to a wave that was, in a 
classical-like, probabilistic sense, fully outside the initial state or the inner region at intermediate 
times, i.e., to a completely reconstructed state, whereas the decay during the exponential regime is 
due instead to a non-reconstructed wave. At short times quantum interference between regenerated 
and non-regenerated paths is responsible for the deviation from the exponential decay. We may 
thus conclude that state reconstruction is a "consistent history" for long time deviations but not 
for short ones. 
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Introduction. The decay from unstable quantum states 
is a recurrent and important tòpic in several fields of 
physics. The exponential decay, which is by far the most 
common type, is surrounded by deviations at short and 
long times . The short-time deviations have been much 
discussed, in particular in connection with the Zeno ef- 
fect, and experimental observations have been performed 
quite recently 0- An exactly solvable analytical model 
for a single resonance that describes these regimes foi- 
lows for a decaying particle of mass m from some general 
considerations of the survival amplitude 0] , 

A(T) = (^ole-^^o) = ^-J^e-* k2T / n A(k)dk, (1) 

(The units employed are ti = 2m = 1.), where the con- 
tour C goes from c + ioo to ic + oo, with c > a small 
constant, along the upper half complex fc-plane and ipo 
is the initial state. If A(k) is given by a single reso- 
nance, A(k) = 2k[D r /(k - k r ) + D*/(k + k*)], and the 
input to evaluate Eq. (J2J) is the value of the complex 
pole k% — e r — iT r /2. It follows from the conditions 
A(T = 0) = 1 and A(k = 0) = that @], 
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[D r w(-y r )+D; W (y_ r )] (2) 



where D r = k r /Ke(k r ), y r = — exp(— Í7r/4)/c r T 1 / 2 , 
y_ r = exp(— Í7r/4)fc*T 1 / 2 and the function w{z) = 
exp(— z 2 )eríc(— iz). The second term in Eq. J2J) doni- 
inates at short and long times compared with the life- 
time. For other analytical or semianalytical examples 
see 0, 0, @ i which show that in general the exponential 
decay may be related to a pole residue and the deviation 
to a linc integral resulting from a contour deformation 
along a steepest descent path in the complex momentum 
plane. The long time deviations may also be associated 
with the boundedness of the spectrum and the short 
time deviations with the existence of the energy moments 
, but clearly a more physical and intuitive understand- 
ing is desirable 0, II] ■ In this respect Fonda and Ghirardi 



in 1972 7] put forward an interesting connection follow- 
ing Ersak's previous insight : " Concluding, we can say 
that the physical mechanism producing deviations from 
the exponential decay law of an unstable system is there- 
fore provided by the regeneration of the initial unstable 
state. In this regeneration process, a very decisive role 
is played by the decayed states, which, as a consequence 
of the evolution of the system, are able to reconstruct 
partially the initial unstable state through a process of 
rescattering" 0- While that paper was centered about 
long time deviations, in fact the arguments given are also 
vàlid for short times, or arbitrary deviations [ïoj . as we 
shall see. Let us decompose the evolution up to the final 
time T into different "paths" with the aid of the following 
complementary projectors 



P = |^o)(M Q = l-P, 



(3) 



at an intermediate time i: 



A(T) = (^ole-^-W^P + Q)e- lHt ' h \i) ) 
= A(T - t)A(t) 

+ {Me- lH{T - t)/h Qe- lHt/h \^). 



(4) 



Except for the last term, A(T) would be exponential 
0, • Since the last term involves a path which goes out 
of the initial state at the intermediate time, the deviation 
from exponential decay is attributed to state reconstruc- 
tion. 

The generality of this argument, however, diminishes 
somehow its explanatory power, and we may rightly won- 
der what is the difference, if any, between the physical 
origin of short and long time deviations, in particular be- 
cause some of the mathematical arguments provided to 
explain the deviations are different for long (Paley and 
Wiener theorem) and short times (behaviour at the ori- 
gin, and existence of energy moments). It also seems rea- 
sonable to demand a quantitative rather than qualitative 
description of the reconstruction process. In addition, a 
key conceptual question should be addressed: Since Eq. 
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(0J is an equation for amplitudes, the language used for 
interpreting it is possibly not fully appropriate, and has 
to be made more precise. As it is well known, Feynman 
paths interfere, so that, in general, regarding a sequcncc 
of events (a history) of the system previous to a final or 
measured event as something that has actually occurred 
is not justified, in the sense that probabilities cannot be 
assigned consistently, in general, to the set of alternativc 
histories. This has been emphasized lately by the consis- 
tent histories interpretation, in which a time sequence of 
events can only be considered a "consistent" (classical- 
like) history under specific conditions, in a nutshell when 
the alternative paths do not interfere [ïïl |. 

By applying the resolution of the identity 1 = P + 
Q at an intermediate time t, the survival probability is 
decomposed, using Eq.(QJ and its complex conjugate, as 

S(T) = \A(T)\ 2 = [PP] t + [PQ] t + [QP] t + [QQ] t , (5) 

where we have introduced the notation, 

[XY] t = (ïPo\e- tH{T - t)/h Xe- mt / h \^ ) 

x (^e^Ye^-W^a), (6) 

(or [XY] for short) in which X and Y can be any of the 
projectors P or Q. It is a simple exercise to see that all 
the terms that contribute to the survival probability in 
Eq. JJJJ can be reduced to products of survival amplitudes 
for different times, 

[PP] t = \A{T-t)A{t)\\ 

[PQ] t = A(T-t)A(t)A*(T)-\A(T~t)A(t)\ 2 , 

[QP] t = A(T-t)*A(t)*A(T)-\A(T-t)A(t)\ 2 , 

[QQ] t = \A{T)\ 2 -A{T-t)A{t)A*{T), 

- A(T - t)*A{t)*A{T) + \A(T - t)A(t)\ 2 . (7) 

The diagonal terms [PP}t and [QQ]t have simple "event" 
interpretations: for [PP]t the particle is initially in |^o), 
it is again in \ipo) at the intermediate time t, and ends 
up in the same state at T; for [QQ]t the particle begins 
and ends in but it is in the complementary subspace 
spanned by Q at the intermediate time t. This is there- 
fore the term associated with initial state reconstruction. 
The nondiagonal terms, however, are quantum interfer- 
ences without classical counterpart. Thus the "histories" 
represented by [PP]t and [QQ]t terms are only "con- 
sistent" if the contribution from the interference terms 
is negligible, i.e., if the survival probability may be at- 
tributed essentially to these diagonal terms only. In other 
words, when the histories are consistent we may plainly 
say that the events have indeed happened in one or the 
other way with certain probabilities, without the need to 
invoke virtual paths and complex amplitudes. There re- 
mains of course an arbitrariness in the definition of "neg- 
ligible" , since the interference terms are often small com- 
pared to the others but rarely zero. One should accept, 



clearly, that the "consistency" or "classicality" of the his- 
tories is not absolute and sharply defined but a gradual 
quality that may however be precisely quantified. 

The models. The exact single resonance model for de- 
cay as the one discussed before, see 0, is possibly the 
simplest one for decomposing the survival since A(t) re- 
duces to a known function. One pays a price though, as 
the single resonance model does not contain an explicit 
form of the Hamiltonian or the initial state. The nu- 
merical calculations presented here correspond to a more 
explicit delta-function Hamiltonian model described be- 
low. Note, however, that all bàsic properties of the sur- 
vival decomposition found for the delta model have also 
been reproduced for the simple single resonance model, 
although the results for the later are not shown to avoid 
unnecessary repetitions. 

In the delta function model the particle is confined to 
the half-axis x > with boundary condition ip(x = 0) 
. and subjected to a delta-function potential, 

V = r)S(x-l). (8) 

In all calculations the initial state is the ground state of 
the infinite well, ip(x,t = 0) = \[2 sin(7rx). The evolved 
wavefunction i\)(x,ï) is calculated by using the resolution 
of the identity based on the eigenfunctions of the station- 
ary Schròdinger equation, delta-normalized in k space. 
Since their overlap with the initial function can be calcu- 
lated in the form of an explicit function, ip(x,t) is reduced 
to a single integral which is performed numerically. The 
calculation may be also carried out by expanding ij)(x,t) 
in terms of the resonant-state basis |^ . 

Consider now a model that represents the effect of a 
detector by adding to the potential in Eq. JHl, an ex- 
ternal absorbing potential — íVq for x > 1 12]. The 
solutions of the stationary Schròdinger equation which 
are right eigenvectors of the Hamiltonian with eigenvalue 
E = q 2 — íVq can be written as 

_J^_ f C x {é^ - < x < 1 

Çq{ ' ~~ (2tt) 1 /2 1 e~ l * x - Se^ x , x>í 1 ' 

where k — (q 2 — iVo) 1 ^ 2 is the wavenumber inside, and q 
the wavenumber outside. For scattering-like solutions, q 
is positive. Note the two branch points of k in the com- 
plex q plane. We shall take the branch cut joining these 
points. Similarly, the root in q = (k 2 + iVç,) 1 ^ 2 is de- 
fined with a branch cut joining the two branch points in 
the k plane. The matching conditions determine the am- 
plitudes Ci and S. The right eigenvectors of a complex 
Hamiltonian have biorthogonal partners 4> q (x), which are 
left eigenvectors of H with the same eigenenergy E . 
Alternatively they may also be regarded as right eigen- 
vectors of with eigenvalue E* . 

Aside from the scattering solutions, there may be lo- 
calized solutions for complex qj and positive imaginary 
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FIG. 1: ln of the survival probability S(t) for different delta 
function strength factors. 



FIG. 3: Decomposition of the survival probability for a final 
time T — 3 where exponential decay holds. n = 5. 
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FIG. 2: Survival probability S(T) and its decomposition as 
a function of the intermediate time t for a final time T — 10 
in the long-time, non-exponential decay regime. r\ — 5. 

part, normalized with respect to their biorthogonal part- 
ners as (uj\ui) = íy. The time evolution in this case is 
obtained from the biorthogonal rcsolution of the identity, 
1 = Ej \uj){uj \ + r o °° dq\(t) q ){^ 3 \. 

Results and discussion. Typical survival curves for dif- 
ferent r\ factors are shown in Fig. ^ (Vo — for the 
time being.) Note the clear transition between exponen- 
tial and long time deviations (The short time deviation 
is not seen in the scale of the figure). 

First we decompose the survival at a large time T in 
which the decay is not exponential for different interme- 
diate times t according to Eq. J5J. Figure |21 shows the 
clear dominance of the reconstructed term. The other 
terms are negligible near the middle intermediate time 
point, 2.5 < t < 7.5, so that the reconstruction is a 
consistent history, in fact the only significant one, in this 
range of intermediate times, whereas the interference and 
non-reconstructed terms are not negligible for extreme 
intermediate times. 

No fundamental difference is observed if instead of the 
survival probability S(T), the nonescape probability N = 
/ °° dx \tl>(x, t)\ 2 is decomposed into terms [PP]t, [QQ]t 
and interference contributions as before, but using now 



FIG. 4: Importance of the interference at short times: De- 
composition of the survival probability for T=0.0001. r\ = 5. 

the new projector P = f Q dx \x) (x\ = Y?jLi\Ej){Ej\> 
where \Ej) are the eigenstates of the infinite well, and 
Q = 1 — P. In practice the sum is only taken up to a 
finit e number of terms until convergence is achieved. The 
calculation reduces then to evaluate matrix elements of 
the general form (E n \e~' lHt / n \Ei) . The resulting figure 
for t] = 5 and T = 3, not shown, is almost identical to 
Fig. 

At variance with the long time case, in the region of ex- 
ponential decay the non-reconstructed term [PP]t dom- 
inates, see Fig. |3 the other ones being negligible at all 
intermediate times. 

In the region of short time deviations the [PP] t , non- 
reconstructed term dominates too. Nevertheless, there 
is a deviation from exponential decay that must be ex- 
plained by the difference with the survival S(T) — [PP]t- 
This difference essentially coincides, and it is numerically 
undistinguishable in our calculations, with the interfer- 
ence term [PQ]t + [QP]t for all intermediate times í, see 
Fig. Therefore, initial state reconstruction in this 
regime cannot be regarded as a consistent history. It is 
an interference effect! 

Fonda, Ghirardi and Rimini |lOj pointed out that the 
state reductions associated with measurement interac- 
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FIG. 5: Efïect of a detector, represented by an imaginary 
potential — iVo outside the internal region limited by the delta. 



present purpose.) 

Thcre is in summary a significant difference between 
short and long time deviations from exponential decay 
regarding the role of state reconstruction, which becomes 
a consistent history for long times but not for short times. 
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FIG. 6: Dominance of the non-reconstructed [PP]-term in 
the survival when the detector is on. 77 = 5, Vo = 1. Compare 
with Fig. H 



tions could hinder initial state reconstruction and make 
the decay purely exponential [ïcH ]. They substantiated 
this idea with a measurement model producing state re- 
ductions at random times. Several models have been 
used afterwards by different authors to show the effects 
of a detector in the exponential decay and its deviations 
|l4l 1 1 rJ | . Here we propose a complex imaginary potential 
as a simple way to represent atòmic detection by fluores- 
cence induced by làser (this is justified in 0). Fig. 
shows that the presence of a weak detector pushes the 
transition between the exponential and long time devia- 
tion to larger times with increasing Vo, showing the grad- 
ual suppression of initial state reconstruction. This can 
also be seen explicitly in Fig. [ïï] where the dominant term 
in the decomposition is [PP] instead of [QQ] as in Fig. 
□ For the same vàlues of Vq, however, the short time de- 
viations are not affected because there is no enough time 
for absorption to occur. Note that the slope of the long 
time deviation increases with Vq. In fact around Vq ~ T 
it actually disappears.(A more extensive analysis of this 
measurement model will be presented elsewhere, in par- 
ticular for strong interactions, which produce important 
changes in the life time, but are not so relevant for our 
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